Conductance of a quantum point contact in the presence of spin-orbit interaction 



(N 
O 
O 
(N 

< 

(N 



Shi-Liang Zhu^^^ ^ 2. D. Wang^-^ *, and Lian Hu^ 
^Department of Physics, University of Hong Kong, Pokfulam Road, Hong Kong, China 
'^Department of Physics, South China Normal University, Guangzhou 510631, China 
^ Department of Material Science and Engineering, University of Science and Technology of China, Hefei, China 

A recursive Green's function technique is developed to calculate the spin-dependent conductance 
in mesoscopic structures. Using this technique, we study the spin-dependent electronic transport of 
quantum point contacts in the presence of the Rashba spin-orbit interaction. We observed that some 
oscillations in the 'quantized' conductance are induced by the spin-orbit interaction, and indicated 
that the oscillations may stem from the spin-orbit coupling associated multiple reflections. It is also 
indicated that the 0.7 structure of the conductance observed in mesoscopic experiments would not 
stem from the spin-orbit interaction. 
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I. INTRODUCTION 

The interaction between the electron spin and its or- 
bital motion, commonly refered to as the spin-orbit (SO) 
interaction or SO coupling, has been known for a long 
time. Recently the effect of SO-interaction on meso- 
scopic transport phenomena and ±he quantum Hall effect 
has attracted growing interest .ErEl Although the interac- 
tion magnitude is small compared to the Fermi energy, 
it may have significant impact on electronic transport, 
particularly in mesoscopic systems where quantum inter- 
ference is extremely important. Using the transfer matrix 
method, Meir et alEl showed that the SO-interaction in 
one-dimensional (ID) non-interacting disordered rings in- 
duces an effective spin-dependent magnetic flux, and then 
any spin-dependent transport quantity can be expressed 
in terms of the same quantity in the absence of the SO 
scattering but with an effective magnetic flux. The adi- 
abatic Berry phase induced by the SO-interaction and 
its effect on the electronic transport were studied fixten- 
sively by Loss et alp Aronov and Lyanda-Geller.0 Per- 
sistent currents in mesoscopic rings induced by the SO- 
interaction was addressed pin Ref.4. Promisingly, as ob- 
served by Morpurgo et-aLjj the geometric phase induced 
by the SO-interactionETEl could induce the splitting of 
the main peak in the ensemble average Fourier spectrum. 
Also interestingly, lifting of the spin degeneracy by the 
Rashba SO-interactionQ was reported experimentally in 
two-dimensional (2D)_electron systems for different semi- 
conductor structures.llEI Lommer et ai. pointed out that 
the Rashba mechanism becomes dominant in a narrow 
gap semiconductor system. Ej The spin splitting energy 
at the Fermi energy in the absence of magnetic field is 
A — 2akF (with a as a SO coupling constant and kp as 
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the Fermi wave vector), which ispabout 1.0 ~ 5.0 meV in 
typical semiconductor materials Ji3 i— in 

On the other hand,|-the observation of the universallj'EJ 
and the nonuniversalEJ quantizations of the conductance 
in a quasi-lD constrication is also remarkable. In a clean 
quasi-lD constrication, if the mean free path is much 
longer than the effective channel length, the conductance 
is quantized in unit of 2Gq with Go = e^/h at zero 
magnetic field, referred to as the universal quantization, 
where-tJftei factor of 2 aries from the electron spin degen- 
eracy.ll3Jl3 Recently, Thomas and co-workers found that, 
in addition to the above quantized conductance plateau, 
there is also a structure at 0.7(200). This so-called 0.7 
structure (nonuniversal or fractional conductance quan- 
tization) appears to be related to spontaneous lifting of 
spin degeneracy in the ID constrication, but the origin 
of it remains as an open question. O Since the 0.7 struc- 
ture may be understood as a zero-field spin splitting with 
an estimated energy of A ~ ImeVfj which is the same 
order as the spin splitting energy induced by the SO- 
interaction, it is nature to ask an important question: 
Is the spin polarization induced by the SO coupling re- 
sponsible for the 0.7 structure? To answer this question, 
we study the effect of SO-interaction on the electronic 
transport through quantum point contacts (QPC) in this 
paper. We find that the 0.7 structure is unlikely to stem 
from the SO-interaction. However, some interesting os- 
cillations in the 'quantized' conductance may be induced 
by the SO-interaction, and may be experimentally ob- 
servable. Moreover, the spin-dependent transport prop- 
erties are of current interest in both fundamental physics 
and applied spin electronics. It is quit intriguing to de- 
velop a method to calculate the spin-dependent transport 
in mesoscopic systems. Also in this paper, genesalizing 
a method established for spin-independent casesE3EJ we 
have done such a job. It isjvorth pointing out that an- 
other theoretical approachllJ may be generalized to cal- 
culate the spin-dependent conductance. 

The paper is organized as follows. In Sec. II A, we de- 
rive the 2D tight-binding Hamiltonian with the Rashba 
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so coupling. In Sec. II B, the recursive Green's function 
technique is presented to calculate the spin-dependent 
transmission coefficient. In Sec. Ill the conductance is 
calculated numerically in a model QPC in the presence 
of SO coupling and/or magnetic field. The finite temper- 
ature effect is also studied. This paper ends with a brief 
summary. 

II. MODEL AND CALCULATION METHOD 

A. 2D tight-binding Hamiltonian with the Rashba 
SO coupling 



H = '^£nm\nm){nm\ - '^{Vnm.n-im\nm){n - 1to| 

nm nm 

+ Vnm,n,n-i\nm){nm - 1| + H.C.), (2) 

where \nm) is a two-component orthonormal set in 
the lattice sites {n,m); = {Unm. + 4i)ro — ijr^ 

with Unm = U{x = na,y = ma), t = TT?'/2m*a?, 
7 = gfiBBz/2; Vnm,n-im = Vxexp(-ieAxa /K) and 
Vnm,nm-i = Vyexp[~ieAya /K) with = Itq - aVs, 
Vy — tTQ + a'ri and a' ~ a /2a. The quaternion ba- 
sis {t^} is defined by the 2x2 unit matrix tq and 
T; = —iai {I = 1,2,3) with ct; as the Pauli matrices. 
Ax or Ay is evaluated at the middle point between sites 
(n, m) and (n -|- 1, m) or (n, m + 1). 
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FIG. 1. Schematic diagram of the system, 
injected from the left lead. 



Electrons are 



In a magnetic field, the 2D Hamiltonian for electrons 
(of effective mass-jm* and charge e) with the Rashba SO- 
intcraction readsLi 



H 



2m* 



a 



+ -((T^rij^ 



U{x, y) 



B, 



(1) 



where 11 = p — |A is the canonical momentum, U{x, y) 
represents the spin-independent potential, A is the elec- 
tromagnetic gauge potential with B = V x A relating it 

to the magnetic fields, and g, /i_B,cr are the g factor, Bohr 
magneton, Pauli matrices respectively, a is the Rashba 
SO-interaction constant which is determined from an ef- 
fective electric field along z direction given by the form 
of the confining potential in the absence of an inversion 
center. Now we choose a discrete square lattice, on which 
points are located at x = na and y = ma, with n, m as 
integers and a as the lattice constant. In terms of the 
quaternion, the one-electron tight-binding Hamiltonian 
with the Rashba SO coupling can now be parameterized 
as 



B. Model and the recursive Green's function 
technique 

To describe spin-dependent transport properties 
through a specific mesoscopic structure, we consider a 
2D structure composed of three different regions ( Fig.l). 
The shadowed central zone is a mesoscopic structure 
where the SO coupling and/or a homogeneous perpen- 
dicular magnetic field are present. To simplify scatter- 
ing boundary conditions, we assume that the structure 
is connected to both sides with semi-infinite ideal leads 
without SO-interaction and magnetic field. In our lattice 
model atoms are at the sites of a square lattice (na, ma) 
( with a to be set equal to unit, m = 1, 2, • • • , M .) in 
the shadowed zone and two ideal leads. As a result 
the system is described by a one-electron tight-binding 
Hamiltonian similar to Eq.(^. For a homogeneous mag- 
netic field B = (0, 0, B) perpendicular to the 2D-plane 
(A = {—yB, 0, 0) in the Landau gauge), we can write 
the Hamiltonian for the shadowed central zone as 



N 



n=l 



N-l 

E 

n=l 



Hc^Y. W)iin{n\ +Y,[\^)Hn,n+i{n+l\ + H.C], (3) 



where \n) is the set of 2M kct vectors belonging to the 
nth cell, and 



/em -Vy+ 
-Vy e„2 -Vy^ 
-Vy e„3 



(4) 



{Hn,n+l)pp' — -V^e 



{p-,v = 1,- 



■,M). 
(5) 



Here /3 — Ba^ with 0o = hc/e as the magnetic flux 
quantum. The forms of Hamiltonian for the two ideal 
leads are the same as Eq.(^), but with different summing 
regions (— oo < rt < 1 at the left lead and N < n < oo 
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at the right lead) and different values of the parameters: 
Vx = Vy = tTo, Enm = ^tTQ, B = 0. Moreover, it seems 
acceptable to assume that the coupling between the ideal 
leads and the mesoscopic structure takes simply the form 
Hq^i — Hn^n+i = —tSpp'.lS 

Following the recursive Green's function tecimiaue. 
developed for some spin-independent systems ,llail3't3 
we may generalize the method to the nontrivial spin- 
dependent cases. The spin-dependent transmission t^'^ 
for the incident channel {I, a) and out-going channel 
(j, u') can be found readily (see the appendix A) . At 
a finite temperature T, the conductance through a 2D 
mesoscopic structure is given by the Landauer-Biittiker 
formulaEJ 



G 



'Go 



E 

lj,<Ta' 



, df{E,T) 



(6) 



where Tj'f 



\tfj' I is determined from Eq.(A3), and 
f{E, T) = [1 + exp{E - Ef)/kBT]-^ is the Fermi-Dirac 
distribution function with fcs as the Boltzmann constant. 
Obviously, Eq (|) reduces to G = GQY.{ij,aa'}^i^]'' ^-t 
zero temperature. 

In the following, we focus on a QPC structure with the 
Biittiker saddle-point potential 



explained by Landauer-Biittiker formula of quantum bal- 
listic transport: for a nonmagnetic QPC, the spin-up and 
spin-down electrons make the same contribution and then 
the unit of the conductance quantization is 2Go. How- 
ever, if the Rashba SO-interaction is present in the QPC, 
the spin degeneracy of conductance electrons would be 
broken. Whether the quantization of the conductance is 
affected by the SO-interaction in a meaningful sense? 
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+ -m*uj'ya'[m-{M + l)l2Y 



(7) 
.0 



which is a practical candidate to describe a real QPC 
uJx,y indicate the strength of the lateral confinement. 
Note that the well-pronounced quantized plateaus occur 

if UJy >U!x- 

Before the end of this section, it is worth pointing out 
that the recursive Green's function method described in 
this section can be straightforwardly generalized to 3D 
mesoscopic structures, and thus is quite useful in study- 
ing the spin-dependent transport of many different meso- 
scopic structures which are of current interest, such as 
the devices with quantum Hall effect, giant magnetore- 
sistance, tunnel magnetoresistence. 



III. NUMERICAL RESULTS 

We now calculate numerically the conductance of the 
QPC in the presence of SO-interaction and/or a magnetic 



field. 

7, TV 



We 
= 15 



2re present the results for the system of M 
\ and t = 1. 



A. The effects of SO interaction 



FIG. 2. The conductance G versus the saddle energy Vb ( in 
unit of El) in the presence of SO-interaction for /3 = ksT — 0, 
Ef — 1.21, LOy = 3lux- The inset shows an enlargement of the 
curve around the arrow. 

Figure 2 shows the conductance G as a function of 
Vb for the SO parameters a' — 0, 0.005, 0.01, and 
0.015. Other parameters are Ef — 1.21, Uy — 3ujx- 
Two effects of the SO-interaction on the conductance are 
seen, (i) The SO-interaction induces some oscillations 
on the plateaus, and the stronger the SO-interaction is 
the larger the oscillation amplitude is. Remarkably, some 
sharp dips exist on the oscillation for the strong SO cou- 
pling, for example, the height of dip is almost 2Go for 
a' ~ 0.015. (ii) A tiny plateau exists near 0.7(2Go) for 
a weak SO interaction a ~ 0.005 ( We will address the 
temperature effect on this tine plateau later). But it be- 
comes a sharp dip for the strong SO-interaction (see the 
inset of Fig. 2). The conductance oscillations as well as 
the sharp dips in the plateaus may be observable in future 
mesoscopic experiments. 

We now attempt to understand the above numerical re- 
sults heuristically. The conductance oscillations appear 
to stem from the multiple reflections with the SO cou- 
pling. By taking into account the multiple reflections, 
the conductance of the QPC may be rewritten asE3 



It is well known that the conductance through QPC 
is quantized with the unit of 2Go as a function of the 
saddle energy Vq. The conductance quantization is well 



G — Go 



^ 1 - 2Rj„cos{2kjaLx) 



i?2 ' 



(8) 
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where the reflection probabihty Rja^ for channel {j, a) is 
given by 



Rja- — 



(9) 



with Vj ( as the value of velocity in the leads ( 
scattering zone). = dEj^/hdk is determined by the 
dispersion relation for the electron states in the QPC, 
which is approximately written as 



.J h^k^/2m* + (j - l/2)hLUy + Vo (10) 

(independent on a) in the absence of the SO coupling,Ei 
but as 



n\^j2m* + aa\kj„\ + (j - l/2)hLjy + Vq (11) 



in the presence of SO coupling. Here we only consider the 
energy spectrum at x = since the transport properties 
for the Biittiker saddle-point potential depend mainly on 
the narrowest part. Note that the quantized plateaus are 
well-pronounced in the absence of the SO coupling (as 



shown in Fig. 2) since v'^ = h kj/m* 



and thus Ri 



for a propagating model in the Biittiker saddle-point po- 
tential is negligible. However, the additional Rj^ may 
be induced from the SO coupling because as seen from 
Eq.(|ri) Vj^ in the presence of SO coupling is different 
from that in the absence of SO coupling. Therefore, the 
oscillations may stem from the SO coupling since Eq (||) 
predicts the transmission oscillations, provided ihat Rj^ 
induced by the SO coupling are not too smalLcj More- 
over, we can understand from Eqs. ^ and ^ that a 
dip in the conductance may appear for a relatively larger 
-RjCT when an electron in the channel (j, a) is in a prop- 

0. It is interesting 



agating state {kj^r > 0) with Vj^ 
to note from Eq.(ll|) that = 



for k^ 



m*a/h?' 



which corresponds the minimum of E"^^ as a function of 
momentum. When Vq varies in the scale of magnitude 
^ akj, Vj^O for the same uj"-channel as in the absence 
of SO coupling. This clearly shows that a larger Rj^ may 
appear in a narrow region in the presence of the SO cou- 
pling, leading to a dip. Our numerical results in Fig. 2 
agree qualitatively with the above arguments. 

We now look into whether the interesting behaviors 
seen in Fig. 2 are sensitive to the set of parameters. 
Firstly, we plot Fig. 3(a) and (b) for two other differ- 
ent ujy/ujx ( the other parameters are the same as those 
in Fig. 2). The plateaus are not well pronounced for 
ujy — l.bcux even if a' = 0(Fig.3(a)). In this case, the 
SO-induced oscillations are not so distinctive with the 
conductance curve for a' = 0, but the sharp dips still 
appear in the conductance. For tOy = 6ujx ( Fig. 3(b)), 
the width of the plateaus is widened. Comparing the re- 
sults in Fig. 2, Fig. 3(a), and (b), we found that the well- 
pronounced plateaus occur for larger ratio ujy/ujx, which 
is due to the Biittiker saddle-point potential. We also 
noticed that the effects induced by the SO coupling are 



similar to those in Fig. 2. Secondly, we plot Fig. 3(c) and 
(d) for two other different Fermi energies (the other pa- 
rameters are the same as those in Fig. 2). It is seen that 
the main effects induced by the SO coupling are still sim- 
ilar to those in Fig. 2. Therefore, we may conclude that 
the essential feature induced by the SO coupling in Fig. 2 
is not quite sensitive to the Fermi energy and the ra- 
tio of LOy/uJx, but is indeed sensitive to the SO coupling 
parameter a' . 
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FIG. 3. The conductance G versus the saddle energy Vo 
in the presence of SO-interaction for j3 = fcsT = 0. (a) 
Ef = 1.21, LOy = l.^LUx\ (b) Ef = 1.21, LUy = Glu^; (c) 
Ef — 2.0, ujy = Zlox\ (d) Ef — 1.0, LOy = Zlox- 

The temperature dependence of the conductance as a 
function of Vq can be seen from Figs. 4(a) and (b). Firstly, 
in certain temperature regions, the oscillation induced by 
the SO-interaction disappears when the temperature in- 
creases, meanwhile the quality of the quantized plateaus 
is improved. However the quantized plateaus are de- 
stroyed when temperature increases further. The mecha- 
nism for the destruction of the quantizedjplateaus by the 
finite temperature is energy averagingJ13 Secondly, the 
tine plateau and the dip near 0.7(2Go) are destroyed by 
energy averaging before the integer plateaus disappear. 
However the 0.7 structure observed by Thomas et al. is 
still observable even at a temperature when all the inte- 
ger quantized plateaus disappear.tJ This obvious differ- 
ent temperature effect on the quantization conductance 
suggests that the 0.7 structure should not stem from the 
SO-interaction. Finally, the quantized plateaus disap- 
pear at /cbT ~ 0.1 for a' = 0.005, while the plateaus 
is still obvious at this temperature for a' — 0.01 (the 
plateau disappears at a higher temperature ksT ^ 0.12). 
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Fig. 4 Zhu et al.. J. Appl. Phys 




FIG. 4. The conductance G versus Vb at finite tempera- 
tures for different SO-interactions: (a) a' — 0.005, with the 
inset showing an enlargement of the curve around the arrow; 
(b) a' — 0.01. Other parameters are Ef = 1.21, Uy = Su^- 

B. The effects of magnetic fields 

We now discuss the effect of a perpendicular magnetic 
field on the conductance quantization: the Zeeman effect 
and the Peierls phase factor. In many theoretical calcula- 
tions, the Zeeman effect is usually ignored for simplicity. 
Here we consider both the Peierls phase and the Zeeman 
effect. Figures. 5(a) and (b) show the field-dependent con- 
ductance as a function of Vq for a' — 0, 0.01, respectively. 
From Fig. 5, one of the characteristic features of QPC's in 
a magnetic field is that conductance steps begin to appear 
at odd integer multiples of Go at /3 > 0.09, due to the 
lifting of the spin degeneracy by the Zeeman effect. An- 
other feature is that the width of the plateaus is widened 
when compared to the (3 = case. On the other hand, 
the number of effective subbands is decreased when the 
magnetic, field is applied, and actually it is proportional 
to 1//3E3 Moreover, figure 5(b) shows that the oscilla- 
tion as well as the sharp dips are destroyed at high fields. 
Thus the quality of the quantization is improved when 
the magnetic field is applied. Certainly, the effect of the 
SO-interaction is suppressed by a high magnetic field. 
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FIG. 5. The conductance G versus Vb in the presence of 
magnetic field, (a) a' = 0, fcsT = 0; (b) a' = 0.01, fcsT = 0. 
Other parameters are Ef = 1.21, ujy = Suix- 

IV. CONCLUSIONS AND DISCUSSIONS 

The SO-interaction is equivalerd. _1x) a momentum- 
dependent effective magnetic field.El'0'l3 However, com- 
paring with a genuine magnetic field, there exist some 
essential differences in the energy spectra of ^ctrons 
induced by an effective field or a genuine field,Ll which 
induce several significant different effects on the conduc- 
tance. Firstly, spin split energy A induced by the SO- 
interaction at |fci = 0| is zero, while the spectrum induced 
by a genuine field is determined by the nonzero Zeeman 
energy gfJ-Bz- This is the reason why the conductance 
steps do not appear at odd integer multiples of Go in 
the presence of SO-interaction but zero magnetic field, 
noting the appearance of this kind of quantized plateaus 
is the essential feature induced by a genuine field. Sec- 
ondly, in the presence of SO coupling, there exist bumps 
( a nonmonotonic portion) in the spectrum as a func- 
tion of transverse momentum, while the spectrum is a 
monotonous function of the momentum in the presence 
of a genuine magnetic field. Therefore, the reflections at 
the scattering region in the presence of SO|-Gpupling may 
induce by the negative velocity for fc > 0Jl3 It provides 
us another understanding of the transmission oscillations. 
Finally, the SO-interaction can remove the spin degener- 
acy in the band structure but could not lead to an over- 
all spin polarization.Q Consequently, the 0.7 structure is 
unlikely induced by the SO-interaction. On the other 
hand, the results addressed in the presence of magnetic 
field may be restricted within weak field since the non- 
interacting electron model used in the recursive Green's 
function is invalid in a very strong magnetic field. 

We are now concerned with the possible experimen- 
tal test of the above SO-induced properties. In unit 
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of t, a' = m*aa/h'^, which is estimated to be (0.54 ~ 
5.4) X 10^'^ for some typical semiconductors [ here we 
choose the parameters of InAs, where a ~ 0.608 nm, 
a ~ (1.0 ^ 10.0) X 10"^° ev cm, andj^* - 0.067 rrie 
with me as the mass of a free electron.t2t] . Thus the os- 
cillations induced by SO coupling may be experimentally 
observable, depending crucially on the value of a' . 

Finally, we wish to make a few remarks on why the 
SO-induced effects on the quantized conductance have 
not been observed in experiments so far. Firstly, the ex- 
perimental observation on the SO-induced phenomenon 
would crucially depend on the value of a', but a' in some 
typical semiconductor systems may be relatively small 
based on the above estimation. Thus it is helpful to in- 
crease the coupling coefficient in experiments. As indi- 
cated by Heida et al|lJ two methods may be used to 
increase the coefhcient a in mesoscopic experiments: one 
is to increase the electron density Hs in QPC since a is 
proportional to Ug ; and the other is to vary the gate bias 
which controls the magnitude of the SO-interaction. On 
the other hand, a larger effective mass of charge carriers 
would also be helpful for observing the SO-induced oscil- 
lations in experiments. Secondly, the SO coupling in the 
leads has been neglected in the present study. This ap- 
proximation may enhance the oscillations observed above 
although the RasUja SO coupling in QPC is larger than 
that in the leads Ea It is extremely difficult to include the 
SO coupling in the leads using the present method be- 
cause the analytical wave function in the leads should be 
given; while it appears to be unsolvable in the leads with 
the SO coupling . Nevertheless, with the advancement of 
nanotechnology, it is possible to fabricate a mesoscopic 
structure where the SO coupling is important only in the 
QPC (i.e., neghgible in the two leads). In such systems, 
it could be easier to observe the SO-induced effects pre- 
dicted here. 

In conclusion, we have generalized a recursive Green's 
function approach to calculate numerically the spin- 
dependent conductance in mesoscopic structures in the 
presence of SO-interaction and magnetic field. The ef- 
fect of SO-interaction on the conductance of a QPC is 
studied in detail. Some interesting oscillations in conduc- 
tance induced by the SO-interaction have been observed 
numerically, which may be tested by future mesoscopic 
experiments. 
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APPENDIX A: 

In the appendix we present the recursive Green's func- 
tion approach to calculate the spin-dependent transmis- 
sion tf^'^ . The electron wave functions in the strip n of 
the leads can be expressed as a linear combination of 
eigenf unctions of a straight infinite lead at a given en- 
ergy E. If the incident electron from the left lead is in 
the (/, a) channel with I and a as the subband and spin 
indices, the scattering states are represented as 



\n)ieft = e''""<l)k,a + '^r/'e ' (Pk,a' 



) right — 2_^tlj e (pkja'. 



(Al) 
(A2) 



where t^''^ (''zj '^ ) transmission (reflection) ampli- 

tudes for the incident channel (Z,cr) and out-going chan- 
nel (j,a'), and the wave functions (/)fc^o- are given by 
(f'k.a = (t>ki ®xi^z) with 



2 nj wjm irjM 

fc, = \l , r . . (s^"^ . .r , -, ^ • • ■ : s^'T' „^ , -, ; ■ ■ ■ , sm- * 



V Af + r M + 1 
and 

' 1 



M + r 



'M + r ' 



= ( ) ' =- \ 1 



Here the sign Tr denotes the transposition of matrix. The 
number of open channels and the value of kj are deter- 
mined by the Fermi energy from the energy dispersion in 
ideal leads 



E = At - 2t{coskj + cos 



■KJ 



M+1 



)■ 



Only if kj is real do we have a propagating state; we have 
an evanescent state when kj is imaginary. All real kj 's are 
positive, because we have already considered separately 
the states propagating to the right (incident and trans- 
mitted wave) or to the left (reflected wave). Following 
the recursive Green's function technique, we flnd 



where 

_ 1 dE 

~ hdk'j' 

M 

1=1 



with 



2 Inp lirp' 

sin sin 

M + 1 M+1 M+1 



(p,p' = l,---,M), 



6 



C(+) 



1 





1 



The set of vectors (pj^.^, are the duals of the set (t>k,a, 
defined by 4>^.„,4'kicr — Sji6a'a- Gn+i,o is the retarded ^ 
Green's function for the scattering region between two 
ideal leads, which can be obtained by a set of recursion ^ 
formulas in a matrix form,cil 

G„'+i,o - .g"'+^i/„'+i,„'G„',o, (0 < n' < N) (A4) 

g"'+i = [£;-i?„,+i-A"']-\ (A5) 
(A"')pp, = e^'^^P-P'^''V+{g"')pp,V,, = 1, • • • , M) (A6) 

by iteration starting form (7'^ — Go.o = (E — Hq)^^, where 
Ho = Ho~ tF{<j,r,), Hi = Hi {1 < I < N), Hn+i = 
Hn^i — tFijTout) with [a out) as the spin state of the 
incident (outgoing) electrons, and 1 



M 



(A7) 



It is worth pointing out that each element of Green's 
function in the above recursion formula is given by a 
quaternion number iru-tJie presence of the SO-interaction 
or the Zeeman effectES In the numerical calculations of 
Gat+i^o and 6((t) in Eg. ([A^) , we have used explicitly the 
multiplication table for quaternion number. 
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